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A vortex in a trapped Bose-Einstein condensate can experience at least two 
types of instabilities. (1). Macroscopic hydrodynamic motion of the vortex 
core relative to the center of mass of the condensate requires some process to 
dissipate energy. (2). Microscopic small- amplitude normal modes can also 
induce an instability. In one specific example, the vortex core again moves 
relative to the overall center of mass, suggesting that there may be only a 
single physical mechanism. 

PACS numbers: 03.75.Fi, 05.30. Jp, 32.80.P. 



1. INTRODUCTION 

The remarkable recent experimental creation of Bose-Einstein conden- 
sates in trapped low-temperature alkali gasesa'BB has generated great inter- 
est in the possibility of vortex states. Such vortices have been widely studied 
in superfluid 4 He,L§til but no clear experimental evidence demonstrates the 
existence of a vortex in a trapped Bose condensate. Theoretical work has 
concentrated on the critical angular velocity S7 c i for vortex creation,! the 
normal modes of a condensate containing a vortex,QtJ'til't3 and general con- 
siderations of stability.0 

The present work studies two specific types of instabilities of a vortex in 
a trapped Bose condensate. The first (Sec. 2) is a hydrodynamic instability 
involving the macroscopic motion of a vortex relative to the background 
condensate, including the effect of the nonuniform condensate density, which 
has not previously been incorporated. The second instability (Sec. 3) 
arises from the microscopic internal oscillations of the vortex, which has been 
studied in a particular geometry by Dodd et al.[i@ This latter behavior is 
especially clear in the weak-coupling limit, where the normal modes of the 
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vortex differ qualitatively from those of a vortex-free condensate. 



2. HYDRODYNAMIC INSTABILITY OF A VORTEX 

The behavior of a trapped condensate depends crucially on the num- 
ber iV of atoms in the condensate. Consider an axisymmetric trap with a 
potential 

V tx = ^M{u 2 xP 2 +u 2 z z 2 ), (1) 

where M is the atomic mass and (p, (ft, z) are the familiar cylindrical polar 
coordinates. The radial and axial oscillator lengths d± = \/H/Mlo± and 
d z = \JfijMuj z characterize the condensate's dimensions for an ideal trapped 
Bose gas (the corresponding volume is of order d\ = d\d z ). The short-range 
two-body interaction potential may be written as 

V(v)^g5^{v), (2) 

where g ~ 4irah 2 /M relates the coupling strength to the s-wave scattering 
length a (here assumed positive). The nonuniform condensate wave function 
ift(r) obeys the Gross-Pitaevskii (GP) equation0'El 

(T + V tl + V H )iP = ^, (3) 

where T = —h 2 V 2 /2M is the kinetic-energy operator, Vff(r) = gN\ift(r)\ 2 = 
ATraNh 2 \ift(r)\ 2 /M is the Hartree potential energy of one particle with the 
remaining particles, and \x is the chemical potential. 

Typically, the dimensionless ratio a /do is small (of order 10~ 3 ), but the 
relevant dimensionless parameterE-3 Na/do varies linearly with the conden- 
sate number N. If Na/do is small, then the interactions are weak, and the 
condensate acts like a nearly ideal Bose gas with characteristic radial and 
axial dimensions d± and d z . In the opposite limit Na/do ^ 1 5 the repul- 
sive interactions predominate, and the condensate expands well beyond the 
harmonic-oscillator lengths. In this Thomas-Fermi (TF) limit, the kinetic 
energy is negligible, and the condensate density follows from the GP equation 

N\^\ 2 = g- 1 (n - V tr ) = n(0) - |J- - ^ , (4) 

where n(0) = fi/g is the central density and R± and R z are the radial and 
axial condensate dimensions with R 2 a /d 2 a = 2n/huj a S> 1 for a =_L and 
z. In the TF limit, the normalization condition Jdy|^| 2 = 1 yields the 
dimensionless parameter 
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where Rq = R\R Z characterizes the TF condensate volume. The repul- 
sive interactions also introduce yet another important length [the coherence 



length £ = 1/ y/8irn(0)a], which here determines the vortex-core radius; in 

'o - 



the TF limit, it obeys the relation £Rq = ^o> implying the set of inequalities 



£ <C d < -Ro- 
lf the trap rotates with angular velocity Q, the relevant "free energy" 
is F = E — £IL Z , where E is the energy of the condensate and L z is its an- 
gular momentum. The phase S of the condensate wave function determines 
the velocity v = (h/M)VS, and the GP equation (|3|) provides an explicit 
expression for the TF free energy 

F^jdV (±Mnv 2 + nV tI + \gn 2 - MnQz • r x v) , (6) 

where spatial derivatives of the density are neglected. 

This integral for F provides a variational expression in the rotating 
frame, and the classical velocity potential $ provides a convenient approx- 
imation for the phase S. The simplest case is a uniform fluid with density 
n per unit length in a long circular cylinder with radius R. When a vortex 
with circulation k = h/M is added to the system at a distance xqR < R 
from the symmetry axis, the change in the free energy AF isc 



AF 



Mk 2 u 
4tt 



]n(f) + l^l-*8)-£(l-*§) 



(7) 



where O = k/2ttR 2 = fi/MR 2 is a characteristic angular speed. In the 
absence of dissipation, the vortex executes a circular orbit at fixed radius 
under the influence of its opposite image at a distance R/xq > R from the 
axis. In the presence of dissipation, however, the vortex moves to reduce its 
free energy. For f2 < Qq, the free energy AF decreases monotonically with 
increasing xq, so the vortex simply spirals outward and annihilates with its 
image. In contrast, if > Qq, the free energy near the axis increases with 
increasing xq, so that a vortex near the axis tends to return to the center of 
the container. This situation is merely metastable for £1 < O c i = fio ln(i?/£), 
since the free energy at the center is higher than that at the wall (with 
a barrier at some intermediate distance), but the vortex becomes a true 
equilibrium state when > Q&. 

The TF radial profile density n = n(0)(l — p 2 /R 2 ) provides a more re- 
alistic description of a nonuniform rotating fluid in a long circular container. 
The mean density per unit length n is half the central density n(0), and the 
corresponding vortex-induced change in the free energy is 



AF = — {l-x 2 ) 



■ (8) 
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An expansion for small x$ shows that a central vortex is stable for Q, > 
O c i = n [21n(i2/O-l], unstable for < n m = O [ln(i?/^) + \] = |ft c i + §, 
and metastable for $l m < Q, < fi c i. Thus, inclusion of a realistic density 
profile affects the detailed form of the free energy (the metastable region is 
considerably reduced), but the qualitative picture remains unchanged. 



3. STABLE AND UNSTABLE MODES OF A VORTEX 

Rokhsar0 has argued that a vortex in a nearly ideal trapped condensate 
is unstable because the relatively large core (of order d±) implies a bound 
state in the core. In particular, he noted that Dodd et alJa found an un- 
stable solution of the Bogoliubov equations for a condensate containing a 
singly quantized vortex (this unstable mode is distinct from the rigid dipole 
modes that necessarily oscillate at the trap frequency oj± independent of the 
interaction parameter Na/do). 

To clarify the physical meaning of this instability, it is helpful to study 
the relevant solutions of the Bogoliubov equations 

(T + V tI - \i + 2gN\iP\ 2 ) Uj - gN{i>)\ = hujUf, (9) 
-gN^fuj + (T + V tI - fi + 2gN\^\ 2 ) V] = -hu jVj , (10) 

where Uj and Vj are the normal-mode amplitudes and LOj is the corresponding 
frequency. A state j with positive norm [J dV(\uj\ 2 — \vj\ 2 ) = 1] is potentially 
unstable if Uj < since the creation of a quasipaxlicle in the jth mode lowers 
the energy relative to that of the condensateJi3 For any solution of the 
Bogoliubov equations, the perturbations in the particle density and velocity 
potential are 



n;- = (ru, - ^i) e-***, *$ = + M e-***. (11) 

For a nearly ideal gas, the terms proportional to g./V in Eqs. (^), (^), and 
(|To| ) can be treated in perturbation theory. If the condensate has no vortex, 
then the leading approximation to the condensate wave function is ip ~ xoo 5 
where I ignore the z-dependent part of the wave function and Xn + ,n-{p, 4>) 
is a normalized two-dimensional oscillator wave function containing n+ and 
n_ right mid left circular quanta created by the raising operators = 
a\ ± iajjX^ The lowest excited solutions of the Bogoliubov equations for 
this ground-state vortex-free condensate are the rigid dipole modes with 
frequency u>± = oj± for all interaction strengths. Apart from corrections of 
order g, their detailed form 

u+ « Xio, «-~xqi, v± ss 0, withu; ± = u;_ L (12) 
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follows either by direct construction or from the general explicit dipole so- 
lutions of the Bogoliubov equations^ 

u± = a±V> v ± = i (13) 
where i/j is any solution of the GP equation ([|). The corresponding density 
and velocity-potential perturbations follow from Eqs. dTT] ) 

n' ± ^n oP e^e-^*; $' ± « ^ P e^e"^', (14) 

where no = | Xoo 1 2 is the unperturbed condensate density. They represent 
a circular motion of the rigidly displaced condensate in the positive and 
negative sense, respectively. 

The situation is much more interesting for a singly quantized vortex, 
with condensate wave function ip ~ xio ^ e ^ P e P ^ 2 hi the noninteracting 
limit. The general construction in Eq. (|il]|) now yields 



\ = ( V2x2o\ f u - \ = fXn 
WJ V Xoo J ' \v-J \ 



with uj± = uj± (15) 



so that the + state involves a coherent superposition of a particle and a hole, 
differing from that for the — state. Nevertheless, they both oscillate with 
frequency iv±, and the densitv and velocity-potential perturbations [from 
Eqs. (|TT|)1 have the same for 



v 
ml 



„< ± M „, (, _ I) e ^ e — , *' ± - ^ (, ± i) e«* ' (10) 

where n„ ~ |xio| 2 is the condensate density for the vortex state (this expres- 
sion also follows by an expansion of the condensate density for a small rigid 
displacement of both the center of mass and the position of the vortex core). 

The present case of a singly quantized vortex leads to an additional 
anomalous mode with frequency uo a ~ —uj±, reflecting the single-particle 
transition from the vortex state to the (lower) Gaussian ground state Xoo-^ 
To zero order in the interaction parameter aN/d z , the associated Bogoliubov 
amplitudes are u a « cosh 8 xoo and v a « sinh 9 X02 , which are properly nor- 
malized for any value of the parameter 8. To determine the actual value of 9, 
it is necessary to use first-order perturbation theory, yielding the expressions 



V X02 




and Cjj a ttU)j_\-l + —\ -\ . (17) 



Like the + dipole mode of the vortex, this anomalous mode is a coherent 
superposition of a particle and a hole. Note that the frequency here dif- 
fers from the trap frequency (and hence from the rigid dipole-oscillation 
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frequency) for any nonzero interaction strength.! The corresponding density 
and velocity-potential perturbation are 

n > ' „ ( p - I) e -^ e -^at and & ~ _JL_ (p+T) e-^e-^K 

s/2 \ p) 2^2Mi \ p) 

(18) 

A detailed study for this anomalous mode shows that the position of the 
vortex core is shifted by twice that of the condensate's center of mass. In 
addition, the overlap integral Jd 2 pn' a (x db iy) vanishes identically, so this 
mode is not excited by the dipole oscillation of the center of mass. 

It is striking that both types of instability (hydrodynamic in Sec. 2 and 
microscopic in Sec. 3) involve the relative motion of the vortex core and the 
center of mass;0 they may well be two descriptions of the same physics. 
It will be particularly interesting to study the character and role of this 
anomalous mode for increasing coupling strength. 
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